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Abstract 



We construct explicit <5-bracketing covers with minimal cardinality for the set 
system of (anchored) rectangles in the two dimensional unit cube. More precisely, 
the cardinality of these (5-bracketing covers are bounded from above by 5~ 2 + o(<5~ 2 ). 
A lower bound for the cardinality of arbitrary 5-bracketing covers for d-dimensional 
anchored boxes from [M. Gnewuch, Bracketing numbers for axis-parallel boxes and 

t^j- . applications to geometric discrepancy, J. Complexity 24 (2008) 154-172] implies the 

lower bound 5~ 2 + 0(<5 -1 ) in dimension d = 2, showing that our constructed covers 

t**"* ■ are (essentially) optimal. 

We study also other ^-bracketing covers for the set system of rectangles, deduce 
the coefficient of the most significant term 5~ 2 in the asymptotic expansion of their 
cardinality, and compute their cardinality for explicit values of 5. 

1 Introduction 

Entropy numbers are measures of the size of a given class JF of functions or sets and 
they are frequently used in fields like density estimation, empirical processes or machine 
learning. Good bounds for these entropy numbers, in particular the covering or the 
bracketing numbers, can, e.g., be used to prove bounds on the expectations of suprema of 
empirical processes (as, e.g., Dudley's metric entropy bound), concentration of measure 
results for these suprema, or to verify that a class T of functions or sets is a Glivenko- 
Cantelli or Donsker Class, i.e., that the corresponding jF-indexed empirical process G n 
exhibits a certain convergence behavior as n tends to infinity (cf. [H [20J, [23]). 

They are also useful in geometric discrepancy theory, i.e., in the theory of uniform 
distribution. (Different facets of this theory are nicely described in the monographs [21 

THE ELECTRONIC JOURNAL OF COMBINATORICS 15 (2008), #R95 1 



[3j [9j [16l HE].) In geometric discrepancy theory one tries to distribute n points in a way 
to minimize the "discrepancy" between a given (probability) measure and the measure 
induced by the points (each point has mass 1/n) with respect to some class of measurable 
sets C. If one takes, e.g., the class Cd '■— {Yli=ii^y x i) \xi,---,Xd *= [0, 1]} of anchored 
d-dimensional axis-parallel boxes, the Lebesgue measure X d on [0, l] d , and an n-point set 
P C [0, l] d , then the so-called star discrepancy of P 



d^P) ■- 



sup 
cec d 



x d (C)--\Pnc\ 

n 



is a measure of how uniform the points of P are distributed in [0, l] d ; here |PnC| denotes 
the cardinality of the set P n S. If one substitutes the set system Cd by, e.g., the system 
of all d-dimensional axis-parallel boxes IZd := {H^ifo) V*) I x ii 2/i> • • • > ^ ^ e [^> 1]}; one 
gets another measure of uniformity, the so-called extreme discrepancy 



dlo(P) 



sup 

c&n d 



x d (C)--\Pnc\ 

n 



Certain types of discrepancy are intimately related to multivariate numerical integration 
of certain function classes (see, e.g., [H [9], [HI [T6], HH US]); a well-known result in this 
direction is the Koksma-Hlawka inequality which, written as an equality, reads 



sup 



f 1 - 

/ f ( X )dx--Y,f(t 



d*oo(h 



i "n)i 



where B is the unit ball in some particular Sobolev space of functions (see, e.g., [14j). 

Thus for multivariate numerical integration it is desirable to be able to calculate the 
star discrepancy of a given point configuration {t 1; . . . , t n }, to have (useful) bounds on the 
smallest possible discrepancy of any n-point set, and to be able to construct sets satisfying 
such bounds. 

Algorithms approximating the star discrepancy of a given n-point set up to some 
admissible error 5 with the help of bracketing covers have been provided in (2TJ [22] (see 
also the discussion in [H]). The more efficient algorithm from [22] generates (^-bracketing 
covers of Cd (for a rigorous definition see Sect. [2j) and uses those to test the discrepancy of 
a given point set. The last step raises the task of orthogonal range counting. Depending 
whether the orthogonal range counting is done in a naive way or (in small dimensions) 
by employing data structures based on range trees, the running time of the algorithm is 
of order 

0{dn\B s \) or 0((d+(\ogn) d )\B s \ + C d n(\ogn) d ), 

where Bg is the generated (^-bracketing cover and C > 1 some constant. The cost of 
generating the (^-bracketing cover Bs is obviously a lower bound for the running time of 
the algorithm and is of order fi(d|jBj|). Thus the running time of the algorithm from 
depends linear on the size of the generated bracketing covers. 



THE ELECTRONIC JOURNAL OF COMBINATORICS 15 (2008), #R95 



Bounds on the smallest possible star discrepancy with essentially optimal asymptotic 
behavior for fixed dimension d have been known for a long time (see, e.g, [2"ll3"l®l[T6"l[r5]). 
Nevertheless, they are nearly useless for high- dimensional numerical integration, because 
one needs exponentially many sample points in d to reach the asymptotic range. Starting 
with the paper [13] probabilistic approaches have been used to prove bounds for the star, 
the extreme, and other types of discrepancy that are useful for samples of moderate size 
[SI El El El EH (HI [IH H7J- I n particular, these investigations focused on the explicit 
dependence on the number of points n and on the dimension d. (Of course, probabilistic 
approaches had been used in discrepancy theory before [13], see, e.g., [H[2]. But these 
studies had not explored the explicit dependence on the dimension d.) 

Let us describe these results in more detail: We denote the smallest possible star 
discrepancy of any n-point configuration in [0, l] d by 

0»,d)= r inf , dUP) 

PC[0,l] d ;\P\=n 

and the so-called inverse of the star discrepancy by 

n^e, d) = min{n G N | d^n, d) < e} . 
In [13] Heinrich, Novak, Wasilkowski, and Wozniakowski proved the bounds 

d* 00 (n,d)<cJ^ and n^e , d) < \ C 2 de' 2 ] , (1) 

where C is a universal constant. The proof uses a theorem of Talagrand on empirical 
processes [20l Thm. 6.6] combined with an upper bound of Haussler on so-called covering 
numbers of Vapnik-Cervonenkis classes [12]. (Since the theorem of Talagrand holds not 
only under a condition on the covering number of the set system S under consideration, 
but also under the alternative condition that the ^-bracketing number of S is bounded 
from above by (C5~ 1 ) d , C some constant [2"0] Thm. 1.1], one can reprove ([T]) by using 
the bracketing result [HI Thm. 1.15] instead of the result of Haussler.) 

An advantage of (PQ) is that the dependence of the inverse of the discrepancy on d is 
optimal. This was verified in [13] by a lower bound for the inverse, which was improved by 
Hinrichs [15] to n^((i, e) > code" 1 . A disadvantage of (CD) is that so far no good estimate 
for the constant C has been published. 

An alternative approach via using bracketing covers and large deviation inequalities 
of Chernov- Hoeffding type leads to slightly worse bounds with explicitly given small con- 
stants [3 El [71 El [HI [13] . Let Nn(Cd,5) denote the bracketing number, i.e., the cardinality 
of a minimal ^-bracketing cover of Cd- Then 



<oM)< 



4(lniY[](C d ,£/2) + ln2) 



s 2 



(2) 



see [H Proof of Thm. 3.2]. Thus improved bounds of the bracketing entropy In Nn(Cd, S) 
would lead directly to improved bounds on the inverse of the star discrepancy and of the 
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star disprepancy as well (although its dependence on the entropy cannot be expressed by 
an explicit formula like ([2D, since the corresponding parameter 8 should be chosen to be 
of the order of the star discrepancy; see again [HJ Proof of Thm. 3.2]). 

Attempts have been made to provide deterministic algorithms constructing point sets 
whose star discrepancy satisfies the probabilistic bounds resulting from this alternative 
approach [SI [TJ, [8] . The running times of the algorithms depend on the cardinality of 
suitable ^-bracketing covers; smaller covers would reduce the running times. 

These examples show that for discrepancy theory and its application to multivariate 
numerical integration it is of interest to be able to construct minimal bracketing covers. 

In [HI Thm. 2.7] we derived for fixed dimension d the upper bound 

N u (C d ,8)<p- d + 0(S- d+1 ) (3) 

for the bracketing number of the set system C d . In [11] the bounds 

6- d (l - c d 8) < N {] (C d , 5) < 2^(5 + l)~ d , (4) 

where c d depends only on the dimension d, where proved. Obviously there is a gap between 
the upper bounds and the lower bound. In this paper we prove that in dimension d = 2 
the lower bound is sharp. More precisely, we construct explicit 5-bracketing covers TZs 
whose cardinality is bounded from above by 5~ 2 + o(5~ 2 ); thus 1 is the correct coefficient 
in front of the most significant term in the expansion of the bracketing number N[](C d , 5) 
with respect to o 7-1 . Furthermore, we discuss other constructions in dimension d = 2 (e.g., 
the cover from [22]) and compare them. We conjecture that the lower bound in (|3j) is 
sharp in the sense that N^(C d , 5) = 5~ d + o d (5~ d ) holds for all d; here o d should emphasize 
that the implicit constants in the o-notation may depend on d. We are convinced that 
this upper bound can be proved constructively by extending the ideas we used to generate 
H§ to higher dimensions. 

2 Preliminaries 

Let d G N and put [d] := {1, . . . , off. For x, y G [0, l] d we write x < y if Xi < y% holds 
for all i G [d]. We write [x,y] := Ylieis [ x ii Vi\ anc ^ use corresponding notation for open 
and half-open intervals. We put V x := A d ([0,a;]) and V x>y := X d ([x,y]), where A d is the ef- 
dimensional Lebesgue measure. Similarly, we put Va '■= A d (A) for any measurable subsets 
A of [0, l] d . In this paper we consider the classes 

C d = {[0,x)\xe[0,l] d } and TZ d = {[x, y) \ x, y G [0, l] d } 

of subsets of [0, l] d . The elements of C d are called anchored (axis-parallel) boxes or simply 
corners. The elements of lZ d are called unanchored (axis-parallel) boxes. (Here the word 
"unanchored" is of course meant in the sense of "not necessarily anchored" . ) 
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Let T G {C d , 7t d }. For a given 5 G (0, 1] and A, B G T with A C 5 we call the set 

[A,%:={CG J|iCCCB} 
a 5-bracket of JF if its weight W([A, B]) defined by 

^([A,s]):=y B -Vk 

does not exceed 5. A 5-bracketing cover of JF is a set of (^-brackets whose union is J- . By 
A^[](jF, 5) we denote the bracketing number of J 7 , i.e., the smallest number of 5-brackets 
whose union is T . The quantity In JV[ i (JF, <5) is called the bracketing entropy of .F. In [TT] 
we showed in particular that 

N[](C d ,S) <N {] (R d ,6) < (N(C d ,6/2)) 2 . (5) 

The second inequality was verified by using arbitrary 5/2-bracketing covers of C d of cardi- 
nality A to construct 5-bracketing covers of TZ d of cardinality at most A 2 (cf. [TIj Lemma 
1.18]); that is why we can restrict ourselves to the construction of bracketing covers of C d . 

Let us identify the boxes [0, x) in C d with their right upper corners x G [0, l] d . Ac- 
cording to this convention, we identify the bracket [[0,x), [0,?/)]c d with the <i-dimensional 
box [x,y]. 

If we are interested in (5-bracketing covers of C d with small cardinality it is clear that 
we should try to maximize the volume of the 5-brackets used. The following lemma states 
how 5-brackets of C d with maximum volume look like. 

Lemma 2.1. Letd>2,5e (0, 1), and let z G [0, l] d with V z > 5. Put 



x = x(z,5) := I 1 

Then [x, z] is the uniquely determined 5-bracket having maximum volume of all 5-brackets 
of C d that contain z. Its volume is 

s ^ 1/d " " 



„„=|l_^l__j I V, 

(In the case where V z < S it is easy to see that z is always contained in some (5-bracket 
[0, £) with maximum volume V$ = 5.) For a proof of the lemma see [HI Lemma 1.1]. 
Now we state a "scaling lemma" which we shall use frequently throughout the paper. 

Lemma 2.2. Let 5 G (0, 1) and X = (Ai, . . . , X d ) E (0, oo) d . Let 

$(A) : R d - R d , On, . . . , x d ) ^ (X lXl , ..., X d x d ). 

Furthermore, let S C [0, l] d such that $(A)5' C [0, l] d . Then the smallest number of 
S-brackets whose union covers S is the smallest number of ( (IIi=i Xi) 8) -brackets whose 
union covers <$>(X)S. 
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The proof is obvious since scaling a bracket by applying $(A) implies that its weight 
is scaled by the multiplicative factor Yli=i <V 

Let us briefly recapitulate the construction of a (5-bracketing cover Qg from [8 J in which 
the 5-brackets are the cells in a non-equidistant grid. We do so for two reasons: We want to 
compare the cardinality of Qg with the (more sophisticated) bracketing covers we present 
later, and, what is more important, the construction of Qg can be viewed as a "building 
block" of all these bracketing covers. 

We construct the non-equidistant grid 

r 9 = ({ Xo ,x 1 ,...,x Km }u{o}) d , (6) 

where Xo,x%, ...,x K ($,d) * s a decreasing sequence in (0,1]. We calculate this sequence re- 
cursively in the following way: Put xq ■— 1 and X\ : = (1 — 5) l ^ d . If Xi > 5, then define 
x i+ i : = (xj — 5)x\~ d . If x i+ i < 5, then put k(S, d) := i + 1, otherwise proceed by calculating 

Xi+2- 

Since Qg consists of the cells of Tg, i.e., of all closed <i-dimensional boxes B whose 
intersection with Tg consists exactly of the 2 d corners of B, we have 

\G s \ = (K(5,d) + l) d . (7) 

It was shown in [Sj, that Qg is a bracketing cover (without explicitly using this notion) 
and that 

Furthermore, it was shown that the inequality k(S, d) < \-jzi g\ no ^ s ; an< i that the 
quotient of the left and the right hand side of this inequality converges to 1 as 5 approaches 
0. But to make proofs shorter in what follows, it is better to use the more precise estimate 

n(5,d) = ln(<i)(5 _1 + 0(1) as 5 approaches zero. (9) 

a — I 

It follows directly from the following identities which are easy to check: 

(ln(l-<5))- 1 = -r 1 + 0(l) (10) 

and 

ln(l - (1 - S) l/d ) = ln(<J) - ln(d) + O(S) (11) 

as 5 tends to zero. 

Let us now confine ourselves to dimension d = 2 and use the shorthand k(8) for 
k(S,2). Put Oi(S) := (1 - i5) l/2 for i = 0,1,..., [5" 1 ] - 1. Then in fact, k(S) + 1 is 
the minimal number of ^-brackets of heights 1 — ai(6) whose union covers the stripe 
[(0, ax(S)), (1, 1)]; the 5-brackets covering the stripe are the rectangles [(xi, ai(S)), (x , 1)], 
[(x 2 , ai(S)), (xi, 1)],..., [(0, ai(S)), (x k{S) , 1)]. 

Let us more generally define u(5, t) to be the minimal number of 5-brackets of heights 
1 — ax(S) whose union covers the stripe [(t, ai(S)), (1, 1)] for some t G [0, 1]. We calculate 
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again Xq,Xi, ... as above and determine u(S,i) such that x u ($ t t)-\ > t and [(xi, cii (<$)), (x , 1)], 
[(x2, ai(S)), (x\, 1)],..., [(£, ai(S)), (x w (<s,t)-ij 1)] are ^-brackets whose union covers the stripe 
[(£, ai(S)), (1, 1)]. From the construction of the Xi we see that 

«, = (i-f)^-f(i-f)-^ 1 1 :j 1 1 :g^ 

and that 2j +i < £ is satisfied if and only if 

. In (1 - (1 - 5)^) - In (t(l - (1 - 5)-^) + gl - ^) 

%+ ~ ln(l - 5) 

Thus 

In (1 - (1 - 5) 1 / 2 ) - In (t(l - (1 - 5)- 1 / 2 ) + 5(1 - 5)- 1 ' 2 ) 



u(6,t) 



2 



ln(l - 5) 



(12) 



Observe that for t = we have indeed a; (5, 0) = k(8) + 1. We shall use the numbers 
u)(S, t) for different S and £ to show that the last bracketing cover we present in this paper 
exhibits the (asymptotically) optimal cardinality. 

In the following three sections we present ^-bracketing covers with reasonably smaller 
cardinality than Qg. 

3 The Construction of Thiemard 

Before stating the algorithm of Thiemard to construct a ^-bracketing cover Tg, we want to 
explain its main idea in dimension d = 2. (In [22J the algorithm is discussed for arbitrary 
d.) 

It covers [0, l] 2 successively with 5-brackets by decomposing all rectangles P with 
weight W(P) > 5 into smaller rectangles starting with the rectangle [0, l] 2 . More precisely, 
if P is of the form P = [a, 13] for some a = a p , /3 = (3 P G [0, l] 2 , then it calculates 
parameters 71 = j p , 72 = 7<f satisfying «i < 71 < j3\ and a 2 < 72 < P2 and decomposes 
P into 

Q p = [(ai, a 2 ), (71, fo)] and P p = [(71, a 2 ), (pi, (3 2 )]- 

Afterwards it decomposes P p into 

Q 2 3 = [(7i,a 2 ),(/?i,7 2 )] and P p = [( 7l , T2 ), (ft, ft)], 
resulting in the (almost disjoint) decomposition 

P = QfuQfuP 2 P 

The right choice of 7 = (71,72) ensures W(P P ) = 5 and P p is chosen to become an 
element of the final 5-bracketing cover Tg . 
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The rectangle Q p is of "type 1", the rectangle Q p of "type 2": if the algorithm 

decomposes them, then it chooses 7 X x G (af,7 P ) and j x 2 — j p implying that Q p will 
be decomposed into three, but Q P only into two non-trivial rectangles. 

That is why in the algorithm a rectangle P is described by the triple (P,i,W(P)), 
where i G {1,2} denotes the type of the rectangle. 

Denoted in pseudo-code, the algorithm looks as follows: 

Algorithm THIEMARD 

Input: 5 G (0,1). 

Output: A (^-bracketing cover %. 

Main 

Decompose ([0, l] 2 , 1, 1) 
Procedure decompose {P,j,v) 
Compute 5 P according to ([TBI 
Compute 7 P according to ( fl4l) 
If 5 p v > 5 

For % from j to 2 

Decompose (Q p ,i,S p v) 
Else 

For % from j to 2 

r, : = r, u {gf } 

r,:=r 5 u{[ 7 p ,/5 p ]} 

For each triple (P,j,v) we calculate <5 P G (0, 1) and 7 P G [0, l] 2 as follows: 

73M^J ,fJ=li d = ^r^r ,fJ=2 ' (13) 

and 

7 P = { ttf iH<J ' (14) 

That the resulting set T$ is indeed a 5-bracketing cover was proved in [22]. In Figure 
[T] and [2] one can see the resulting cover Tg for 5 = 0.25 and 5 = 0.05. 

Let us now determine the asymptotic behavior of \T&\ for S tending to zero. In 
Theorem 3.4] Thiemard proved the bound 



\T/j\ < I I , where h 
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ln(l - 5) 
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Figure 1: T 5 for 5 = 0.25. 
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Figure 2: T 5 for 5 = 0.05. 
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This implies \%\ < 2(ln(5~ 1 )) 2 5~ 2 + o(5~ 2 ). We improve this estimate in the following 
Proposition by deducing the correct asymptotic behavior in terms of £ _1 and the exact 
coefficient in front of the most significant term 5~ 2 . 

Proposition 3.1. For a given 5 G (0, 1) we get 

\%\ =21n(2)(T 2 + 0(cr 1 ). 

Proof. From the discussion above (and also from Figure [1] and [2]) we see that Thiemard's 
algorithm decomposes the unit rectangle [0, l] 2 into stripes 

Sf :^[{t i ^Q),{t h l)l i = 0,...,r(5), 

and these stripes again into 5-brackets; here the numbers t, are the x-coordinates of the 
corners of all rectangles of type 1 that appear in the course of the algorithm. More 
precisely, we have to = 1, t r (<5)+i = 0, 

ti+1 ,( 1 -lX i2 ti J( ti -i) 2 ^) ll \ (15) 



tij \\ ^ 

and r(S) is uniquely determined by the relation 

< t T(s) < 5. (16) 

We have 

x 

U-S< t i+1 < U - -; (17) 

both inequalities follow easily from (TT5|) . From (TT6|) and (ITT]) we get 

[r 1 ! - 1 < r(8) < [2r 1 ] - 1. (18) 

Furthermore, we get by simple induction 



<?+i = l-5^t fc , 



fc=0 

which, together with f!16p . results in 

r(8)-l 

5- 1 -5< J2 f k < S ~ X - ( 19 ) 

fe=0 

Let us now calculate the number sjj of (^-brackets of widths t, — tj + i that cover the stripe 
Sg . Since the bracketing problem is symmetric in the x- and y-coordinate, we get from 
the discussion in the previous section 

s$ = k(5) + 1, where k(5) = k(5, 2) as defined in (JS]). 
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(Note that ti — (1 — 5) 1 ^ 2 .) From this we can derive Sg for all i 6 {0, ... , r(5) — 1} via 
"scaling": Lemma l2~2l gives us with the choice A = (t^ 1 , 1) 

s f = K (6/U) + 1 for all i e {0, . . . , r{8) - 1}. 

(Observe that t i+ i/ti = ai(5/ti).) Furthermore, we have trivially sj> = 1. 
Then identity Q) and the inequalities (fT8l) and (fT9l result in 

r(5) r(S)-l t(S)-1 

\T 5 \ = J2sf= J2 (^/U) + l) + l = 2\n(2)5- 1 ]T t. + O^ 1 ) = 21n(2)r 2 + 0(r 1 ). 

i=0 i=0 i=0 

D 



4 Another Construction 

Let us consider another algorithm constructing 5-bracketing covers Zg for anchored rect- 
angles: 

Let again a« = a«((5) = (1 — id) 1 / 2 for i — 0, . . . , C(<5) := |~<5 _1 ] — 1, and a^(s)+i = 0. Put 
Oj := (aj(5), aj(5)) for all i. We first decompose [0, l] 2 into layers 

L^{S):=[0,a i ]\[0,a i+1 ). 

Then, starting with L^°'(8), we will cover each layer L^\S) separately with 5-brackets. To 
this purpose we cover for fixed i G {0, . . . , ((5) — 1} the stripe [(0, aj +i ), (a^, a^)] recursively 
by the following procedure: 

Put Si(5) := {[a i+ i,ai]} and x 1 := a i+1 (5). 
If XjCLi(5) > 0, then define 

x j+ i := max{0, (xja,i(5) - 5)/a i+1 (5)} 

and put 

Si(5) := Si(5) U {[(x j+1 , a i+ x(6)), (xj, a,i(8))]}. 

If XjCLi(5) < 5, then stop the covering procedure. 

It is easy to see that for each i the resulting set Si (5) consists of 5-brackets whose 
union is [(0, Oj+i), (oj, Oj)]. In fact, we see that for i = the Xj, j = 0,1,..., we get 
from the procedure above form the projection of r<5 (defined as in ([6])), i.e., the set 
{x ,X!,. .. ,x K (s),x K (s)+i}, where x K (g) < $ and x K (g)+i =0. Thus \S Q (5)\ = k(S) + 1. Using 
the scaling Lemma l2~2l with A = (o^ , aj 1 ) we deduce that consequently |«Sj(5)| = n(5i)+l, 
where Si := <5/(l — iS). (Observe that a i+ i(8) / cti(8) = ai(Si) for i < ((5).) By symmetry, 
we can cover L^'(j) by 2k(g\) + 1 5-brackets. (Observe that L^(S) is already a 5-bracket.) 
More precisely, using the mapping ref : M 2 — > R 2 , (x,y) i— >• (y,x) we have 

C(S) 
Z 5 = \J(Si(5)Uiei(S l (5))). 

i=0 
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Figure 3: Z 5 for 5 = 0.25. 



The set Z$ is a 5-bracketing cover of [0, l] 2 with 

/C(*)-1 \ 

\Z S \= Yl (2«W + 1) +1- 

In Figure [3] and H] we see £$ for 5 = 0.25 and 0.05. 

From the identity ([20]) . the definition of ((S), and (jHJ) we get 

/CW-i \ 

1^,1 = 2 I £ 21n(2)(Jr 1 + O^- 1 ). 

-i _ x-i 



Now 5j — S — i, hence 



CW-l 



£ V = «*)*-' - «S«^i) = !*-> + 0(0- 



8=0 



Thus 



|Z a | = 21n(2)<J- 2 + 0( ( J- 1 ) 
Altogether, we proved the following proposition. 



(20) 
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Figure 4: Z 5 for 5 = 0.05. 



Proposition 4.1. For 5 G (0, 1) the set of rectangles Z$ constructed above is a 5- 
bracketing cover of C<i . Its cardinality is given by 



'«*)-! 



\Z* 



J2 (2«(«Ji) + 1) +1 = 2 ln(2)r 2 + 0(0, (21) 

i=0 I 

where ((5) = [5 _1 ] — 1, Si — 5/(1 — i5), and K(Si) = K,(Si, 2) as defined in (TJJ). 

5 Re-Orientation of the Brackets 

A positive aspect of the two previous constructions is that (essentially) all brackets in the 
resulting ^-bracketing covers have largest possible weight 5 and overlap only on sets of 
Lebesgue measure zero. But if we look at the brackets in Thiemard's construction which 
have some distance to the upper edge of the unit rectangle [0, l] 2 , then these boxes do 
certainly not satisfy the "maximum area criterion" stated in Lemma I2TT1 The same holds 
for the brackets in Z$ which are close to the x- or the y-axis and away from the main 
diagonal. The idea of our next construction is to generate a bracketing cover similarly as 
in the previous section, but to "re-orientate" the brackets from time to time in the course 
of the algorithm to enlarge the area which is covered by a single bracket. Of course the 
algorithm should still be simple and avoid to much overlap of the generated brackets. 
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Before stating the technical details, we want to present the underlying geometrical 
idea in a simplified way: 

Like the construction in the previous section, our new bracketing cover should be 
symmetric with respect to both coordinate axes. Thus we only have to state explicitly 
how to cover the subset 

H:={(x,y)e[0,l] 2 \x<y} 

of [0, l] 2 . For a certain number p = p(S) we then subdivide H into sectors 



T"':= (l,,)£ff\{(0,0)) 



V £ ^^ U{(0 ' 0)} ' " = 1 -"- 2 "- 



In the same manner as we decomposed in the previous construction the set H into stripes 
[(0, a i+ i), (oj, Oj)], we now decompose T^ into stripes [(0, a i+ i), (a*, Oj)] nT^ 2P \ We do it 
similarly with the sectors T^\ . . . ,T^ 2P ~ 1 \ but we use thicker (and therefore less) stripes 
there. Covering each of these stripes by brackets whose height is exactly the height of the 
corresponding stripe, we see that each bracket has almost the maximum possible area. 
Provided we can avoid to much overlap at the boundaries of the sectors, we thus need 
only a very small number of these brackets to cover [0, l] 2 . 

Let us now state the generating algorithm precisely. We define "discretized" versions 
T^J of the sectors T^ h \ composed of stripes. To this purpose we define for each h E 
{1,...,2*} 

p ih \5) := \h2- p 5- 1 ] - 1. 

(Note that p^ 2V \5) is precisely ((8) as defined in the previous section.) For i — 0, . . . , p( h \8) 
let 

T?\8) : = [(4S(5),o£ ) 1 (5)),(/ i 2^of ) (5),of ) (5))], 

where 

h-i : 



4£m - V I ' - 



-i 



h h 



op \ 1/2 

" 5 



h-l 



and 



«r{S):={l-i^6 



.CO/ 



1/2 



(Here we use the convention to denote for a general function / by / + the function 
/l/-i([o,oo))> where 1a is the characteristic function of a set A. In particular we have 
t^(S) = for all i and a { ^ h){s)+1 (S) = for all h.) We put af\d) := [hl^af \5), of \5)). 
Then 

^ == U *f 'W 

i=0 

can be viewed as a discretized version (discretized with respect to a decomposition into 
stripes) of T"C0. 
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Now for h = 2 P , 2 P - 1, . . . , 1 we cover each stripe T t {h) (5), i = 0, 1, . . . , p {h) {5), of the 
"discretized" sectors T^J by brackets having exactly the height of the stripe T^ (5) in 
the following manner: 

Algorithm RE-ORIENTED BRACKETS 

Input: 5 e (0,1), p e [0, oo). 
Output: A (^-bracketing cover IZg. 

Main 

For h = 2 P to 1 

For i = 0,...,pW(5) 






Ks^KsUiia^Slaf'iS)]} 



For j = l,2,... 

lix 3 >tfl{5) 



x 3+l :=(x 3 a { t\5)-6) + (a { SMY 1 

K 8 := TZ S U {[(x j+l ,a^ x (S)), (x^afV))]} 



Else next i 
lis ■= K 5 U Tei(Tl s ) 

The output set TZs is visualized in Figure [5] and [6] for 5 = 0.05 and 5 = 0.02; there we 
have chosen p = p(5) to be 

ln^" 1 ) - k 



p(5) = max 



(22) 



with k = and c = 1.7. With this choice we get for 5 = 0.25 that p = and consequently 
Z$ = TZs; thus Figure [3] shows TZs for <5 = 0.25. 
Let us now prove the following proposition. 

Proposition 5.1. The output set TZs of the algorithm stated above is a 5-bracketing cover. 
If p — p(6) is a decreasing function on (0, 1) with \ims->op(5) = oo and 2 P = o(5~ 1 ) as 5 
tends to zero, then the bracketing cover TZs satisfies 

\TZs\=5- 2 + o(5- 2 ). 
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Figure 5: K s for 5 = 0.05. 
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Figure 6: K s for 5 = 0.02. 
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Proof. One can check by direct calculation that all rectangles that are added to TZs are 
in fact 5-brackets. The points a\ (5) are lying on the lines - = ^ and the x-coordinates 



dh) 



<h) 



y — 2P 



2P rp{h)_ 



ti+i($) of the left corners of the stripes T^ (5) are chosen in such a way that H C u|! =1 T^ is . 
For given h > 2 and a given i the index 



'h-1. 1 



is uniquely determined by 



. (fc- 1 )/ 



» 



(ft-i). 



arw>cw>«yw. 



and we have 



,""^ = ^LL a <*- ] >, 



O*) 



-r } ^- 



(h). 



» 



Let now u;(<5, /i, z) be the minimal number of (^-brackets of heights a\ (S) — a 4 - +1 (5) that 
we need to cover T> {§). Using the scaling Lemma I2T21 with A = (2 p (ha i (5))~ 1 , (a,i(5))~ l ) 
we see that we have u>(S, h, i) = u{8\ ,t(5, h, i)) as defined in Section [21 where 5^ = ^5 
and, coinciding with the convention from the previous section, 5\ = 5^/(1 — iS^ h '), and 



t(S, h,i) 



h-1 
h 



h-1. 1 
h l ~h 



*-- A^ ,/s 



h-1 



Due to flTJ) we get 



u;(<5, /i, i) 



In (1 - (1 - <f } ) 1/2 ) - In {t(S, h, z)(l - (1 - <f V 1 / 2 ) + <f >(1 - 5f VV2) 



ln(l - £ 



(ft), 



We claim that 



.,:« <) /; . / ) = 2 In ( 1 + i J (<jf >) L + 0(1) as S^ h) tends to zero. 



(23) 



According to ([H]) this is true for h = 1. In general it follows from the inequalities (fTUj) . 
(fTTj) and 



In (^,/ lj i)(l-(l-^>)-V2 ) + (5 W (1 _ (5 W ) -i/2 ) = in A + 1A -l n (2)+l n (^)+0(c5f>). 



We have 



2P /p w (S) 
ft=l \ i=0 



(24) 
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here we have to subtract the last term to avoid double-counting of the (^-brackets on the 
main diagonal of [0, l] 2 . According to (j23p we get 

h=l i=0 V ' ^ ' 

op / 9 \ 

It remains to show that the sum in parentheses is of the form | + o(l) as 5 tends to zero 
(and thus p tends to infinity). But this follows easily from the identity 



mll + I^ =h- 1 -h- 2 JT(-±) k ^ 

fc=0 



h 



k + 2 



□ 



6 Numerical Comparison and Conclusion 

Let us now compare the cardinalities of the different constructions of (^-bracketing covers 
for some values of 5, see the table below. For the construction of TZ$ we have chosen 
P = p(3) exactly as in ( 1221 . Thus 



2 P ^5' 



ln(2) 



5 



-0.4 



o^ 1 ), 



and the conditions of Proposition 15.11 are clearly satisfied. Note that 21n(2) = 1.386294... 
and (2 ln(2)) 2 = 1.921812.... Thus the table underlines the dominance of the leading terms 
in the expansion of the cardinalities of the ^-bracketing covers with respect to <5 _1 . 



5 


0.25 


0.1 


0.05 


0.01 


0.005 


0.001 


0.0005 


0.0001 


\Gs\ 


36 


196 


784 


19321 


77284 


1923769 


7689529 


192182769 


\%\ 


25 


142 


565 


13922 


55575 


1386908 


5546403 


138635574 


\Zs\ 


24 


146 


572 


13962 


55650 


1387292 


5547174 


138639434 


\K S \ 


24 


128 


490 


10888 


42162 


1021122 


4055986 


100514774 


5- 2 


16 


100 


400 


10000 


40000 


1000000 


4000000 


100000000 



Altogether, we provided in this paper an explicit construction of a (5-bracketing cover 
1Z$ of C2 which is optimal in the sense that the coefficient in front of the most significant 
term 5~ 2 in the expansion of \1Z§\ with respect to £ -1 is optimal. 

We compared TZs to its simplified version Zg (which does not "re-orientate" the brack- 
ets) and known bracketing covers from [8] and |22j . 
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We conjecture that extending the idea of construction of 1Z& to arbitrary dimension 
a can generate (^-bracketing covers lZ g of Cd whose cardinality satisfies 



\TZ S '\ — 8 d + Od(5 d ) as 5 approaches zero 

(here Od should emphasize that the implicit constants in the o-notation may depend on 
d), i.e., has the best possible coefficient in front of the most significant term 5~ d in the 
expansion with respect to S -1 . 

We suspect that a rigorous proof of the conjecture might be rather technical and 
tedious. That is why we would find even a rigorous analysis for d = 3 or computational 
experiments for higher dimension quite interesting. 
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